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ABSTRACT 

The X-ray transient MAXI J1659— 152 was discovered by Swift/BAT and it was ini- 
tially identified as a GRB. Soon its Galactic origin and binary nature were established. 
There exists a wealth of multi- wavelength monitoring data for this source, providing a 
great coverage of the full X-ray transition in this candidate black hole binary system. 
We obtained two epochs of European VLBI Network (EVN) electronic- VLBI (e-VLBI) 
CO . and four epochs of Very Long Baseline Array (VLBA) data of MAXI J1659— 152 which 

show evidence for outflow in the early phases. The overall source properties (polar- 
ization, milliarcsecond-scale radio structure, flat radio spectrum) are described well 
with the presence of a compact jet in the system through the transition from the 
hard-intermediate to the soft X-ray spectral state. The apparent dependence of source 
size and the radio core position on the observed flux density (luminosity dependent 
core shift) support this interpretation as well. We see no evidence for major discrete 
53 ! ejecta during the outburst. For the source proper motion we derive 2a upper limits 

of 115 /ias/day in right ascension, and 37 /ias/day in declination, over a time baseline 
of 12 days. These correspond to velocities of 1400 km/s and 440 km/s, respectively, 
assuming a source distance of ^7 kpc. 

Key words: ISM: jets and outflows - X-rays: binaries - stars: individual 
(MAXI J1659-152). 
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1 INTRODUCTION 

Galactic black hole X-ray binaries (BHXRB) are key to un- 
derstanding the power sources of accreting compact objects 
and the physical processes of accretion and jet formation. 
A great number of these systems have been identified as 
transient X-ray s ources, and sometimes they are referred to 
as m icroquasars l|Mirabel et al.lll992l ; iMirabel fc Rodriguez! 
Il994h . owing to their ability to produce collimated, highly 
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relativistic ejecta analogous to those powered by super- 
massive black holes in the centres of active galactic nu- 
clei (AGN). There are two types of jets in BHXRBs: com- 
pact, synchrotron self-absorbed jets that are typical of 
the canonical hard spectral state (e.g. IStirhng et al.ll200ll ; 
iDhawan et al.ll2000l ). and transient ejecta during X-ray state 
transitions from the hard state to the soft state, crossing a 
region in the X-ray luminosity vs . X-ray hardness space o f- 
ten referred to as the "jet-line" dFender et al.ll2004 I2009T ). 
Whether there is a universal jet-line region for all BHXRB, 
what its boundaries are, and most importantly, what the 
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underlying physics is, remain a matter of debat^E Most no- 
tably, we point out that high angular resolution observations 
on milliarcsecond (mas) scales have rarely been available 
in the past, although there has been an increasing num- 
ber of VLBI observations in the past 5 years, thanks in 
part to the more flexible, real-time e -VLBI operations of 
the EVN (see e.g. iTudose et al1l2007l ). Using the e-VLBI 
technique with the EVN to observe transients has the ad- 
vantage of the quick turnaround time, providing more effi- 
cient ways to schedule follow-up EVN and/or VLB A ob- 
servations, as previously demonstrated for various types 
of transients |Miller- Jones et all |2010|; iParagi et al.1 l2010al : 
iMoldon. Ribo fc Paredesll201ll ; lYang et al.ll201ll ). The rapid 
feedback is important for a better understanding of the evo- 
lution of the mas-scale source properties early on during the 
outburst event (especially for transients with no firm iden- 
tification at their discovery). It also gives the opportunity, 
as we demonstrate in this paper, to test observing strate- 
gies and look for closer VLBI calibrators for targets with 
less favourable locations, such as low-declination (< —10°) 
sources for the EVN. The focus of this paper is a new 
BHXRB, which benefited from both. 

On 2010 September 25 Swift/BAT discovered a new 
transient, which was initially identified a s a gamma-ray 
burst and was designated as GRB 100925 A (Ma ngano et al.l 
l2010r >. The MAXI/GSC team detected a hard X-ray tran- 
sient positionally coincident with this GRB candidate, an d 
catalogued it as MAXI J1659-152 (|Negoro et al.l 120101 ). 
Soon it became clear from optical spectroscopy by the 
ESO/VLT X-shooter that the transient showed broad emis- 
sion and absorption lines at zero redshift, indicating that the 
sourc e was likely a Galactic X-ray binary (|de Ugarte Postigq 
2010). The source followed an evolutionary track on the X- 
ray hardness-intensity diagram (HID) and also had vari- 
ability/timing properties typical of BHXRBs, in particu- 
lar type-B and type- C quasi-periodic oscillations (QPO) 
ijKalamkar et alj|201ll ). The X-ray lightcurves from Swift, 
XMM-Newton and RXTE showed irregularly shaped dips 
that recurred with a period of ~2.4 hours, making the source 
the shortest orbital period Galact i c black hole binary candi- 
date known l|Kuulkers et al.|[2012i ; iKennea et al. 2011). The 
two groups that reported these results also showed that the 
companion star is most likely an M5 dwarf, and that the sys- 
tem is likely located at a distance of 7 kpc (but see Sect. 4.1). 
In addition, using the scaling relation between the spectral 
index and the QPO freq uency, the mass of the BH is esti- 
mated to be 20 ± 3 M (|Shaposhnikov et al.1 [2012] ). If this 
estimate is correct, then MAXI J1659— 152 would be the 
most massive stellar black hole known in the Galaxy. 

The apparent high distance from the Galactic plane 
(z = 2.4 ± 1.0 kpc) makes MAXI J1659-1 52 a good can- 
didate for being a "runawa y microquasar" (Yamaok a et al.l 
120121 ; iKuulkers et al.|[2013ll , which received a large kick ve- 
locity during the formation of the black hole in the sys- 
tem. With high resolution VLBI observations one could in 
principle measure the proper moti on even durin g a single 
outburst, as was demonstrated by iMirabel et al.l (|200ll ) in 
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case of XTE J1118+4 80. The claimed mass of 2 ± 3 M 
for MAXI J1659-152 (|Shaposhnikov et alJlioil ) is incom- 
patible with the runaway microquasar scenario, therefore 
constraining the proper motion with VLBI is particularly 
interesting in this case. However, for rapidly evolving out- 
bursts, the short timescales can make such a measurement 
quite challenging. 

Just ~1.5 days after the initial trigger, 
MAXI J1659-152 was detected at the 5 mjy level at 
5 GHz w ith the Westerbork Synthe sis Radio Telescope 
(WSRT) |van der Horst et al.l l2010al ) . This was rapidly 
followed by a series of EVN e-VLBI (|Paragi et al.1 l2010bf ) 
and VLBA target of opportunity observations at 5 GHz. 
This paper summarizes the results of all the VLBI obser- 
vations. We describe the initial e-VLBI experiment with 
the selection of secondary calibrators, and the follow-up 
EVN and VLBA observations in Sect. 2. Data reduction 
and results will be presented in Sect. 3. The compact 
jet model, as a framework for the interpretation of the 
mas-scale structure is presented in Sect. 4. In Sect. 5. we 
compare the mas-scale properties of MAXI J1659— 152 to 
those of other BHXRBs with similar evolution on an X-ray 
hardness-intensity diagram. Conclusions will be drawn in 
Sect. 6. 



2 OBSERVATIONAL STRATEGY AND 

SECONDARY CALIBRATOR SELECTION 

The first observation of MAXI J1659— 152 was carried out 
on 2010 September 30 with the EVN at 5 GHz in re- 
altime e-VLBI mode with the EVN MklV Data Proces- 
sor ("e-EVN"; experiment code RP016A). After the ini- 
tial clock-searching and setup checking, the project started 
at 13:00 UT and lasted 5.5 hours. The participating tele- 
scopes were Effelsberg, Medicina, Onsala, Toruri, the phased 
array Westerbork Synthesis Radio Telescope (WSRT), Jo- 
drell Bank (Mkll), Cambridge and Hartebeesthoek. The 
data rate per telescope was 1024 Mbps, except for Medic- 
ina, Hartebeesthoek (both 896 Mbps) and Cambridge (128 
Mbps, limited by the MERLIN microwave link to Jodrell 
Bank). We used 2-bit sampling and observed both left 
and right circular polarizations. The target was phase- 
referenced to the nearby calibrator J1707— 1415 (2.2 degrees 
away), selected from the VLBA Calibrator listQ The coor- 
dinates used for correlation of the phase reference source 
were a(J2000) = 17:07:20.390556 (±0.25/ cos(5) mas) and 
S (J2000) = -14:15:23. 12829 (±0.5 mas). 

Because of the low declination of the field, the source 
elevations were low, especially for the EVN: the maximum 
elevation for most EVN telescopes was around 20 degrees, 
and it was below 20° at some telescopes for a significant 
fraction of the observing time. This limits the accuracy of 
determining the interferometer phase for the target using 
calibrator measurements. The coherence time is shorter, and 
systematic errors due to small deviations from the correla- 
tor delay model result in errors in the astrometry, as well as 
decreasing the fidelity of the final image. At the frequency of 
our observations, tropospheric and ionospheric phase errors 
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have comparable effect on the data, however at low eleva- 
tions the tropospheric errors will likely dominate. Both at 
the EVN and the VLBA correlators the tropospheric delays 
are determined from a zenith value which is then mapped 
to lower elevations. A small error in the tropospheric zenith 
delay, imperfections in the mapping function itself, or in- 
homogeneities in the troposphere will lead to a significant 
differential phase error between the calibrator and the tar- 
get close to the horizon. This results in correlation ampli- 
tude losses, degraded image fidelity, and poor astrometry. 
The closer the reference source, the smaller the differential 
phase errors between the reference and the target fields. 

Therefore we looked for additional radio sources 
closer to MAXI J1659-152. Unfortunately, this field 
is not covered by the VLA B-array FIRST survey 
l|Becker. White fc Helfandlll995h . which is an excellent re- 
source to loo k for potential se condary calibrators for VLBI 
observations l|Frev et al.ll2008h . Instead, we s elected 8 radio 
sourc es from the VLA D-array NVSS survey l|Condon et al.l 
1 1998T ) that were located within half a degree of the tar- 
get and had total flux densities exceeding 10 mjy. These 
were cross-checked against our 5 GHz WSRT data (from 
ongoing total flux density monitoring). One target, NVSS 
J170003. 28- 145622.9 (hereafter J1700-1456) had a nearly 
flat spectrum. We included this source in the EVN ob- 
serving schedule to verify its compact structure on mil- 
liarcsecond scales. The phase referencing cycle was 90s on 
J1707-1415 and 150s on MAXI J1659-152. Approximately 
every 16 minutes we also observed the secondary calibra- 
tor J1700-1456. Additional calibrators OQ208, J1724-1443, 
J1751+0939 and J1310+3220 were observed for a short 
time to check amplitude calibration of the EVN; 3C286 
was included for calibrating the WSRT synthesis array data 
recorded parallel to the VLBI observations. After the suc- 
cessful detection of both the transient and J1700— 1456 (see 
Sect. 3.) we carried out further VLBI observations, all at 
5 GHz. However, in the case of J1700— 1456 we quickly up- 
dated the coordinates because they were found to be in er- 
ror by about one arcsecond, which is a significant fraction of 
the WSRT phased-array beam at 5 GHz. The second epoch 
was observed with the VLBA (BV070A) on October 2. The 
VLBA data were recorded at the telescopes at a recording 
rate of 512 Mbps, in dual-polarization m ode using 2-bit sam - 
pling, and were correlated with DiFX (|Deller et alj|201lh . 
The third epoch was observed with the e-EVN on October 
4 (RP016B), with a similar setup to the first experiment 
but adding the 40m Yebes telescope to the array. Because 
MAX I J1659— 152 was sti ll detectable with very long base- 
lines (IParagi et al.ll2010bh . additional VLBA observations 
were organized (BV070B-D) on 2010 October 6, 14 and 19, 
until the transient faded below the detection threshold of 
our WSRT monitoring. The time between observing epochs 
was set based on the initial e-EVN results and the X-ray 
spectral and timing properties as observed by RXTE. 
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Figure 1. Naturally weighted 5 GHz VLBA map of J1700-1456 
on 2010 October 2. This source was found to be compact in the 
first e-EVN observations and was used as a secondary phase- 
reference calibrator at the following epochs. The higher resolution 
VLBA observations show a double structure on mas scales, but 
the signal-to-noise ratio was still sufficient for additional phase 
calibration after the residual fringe delay, rate and phase solutions 
were derived using the primary reference source J1707— 1415, 
which lies further away from the target. The peak brightness is 
4.5 mjy/beam. The contour levels are —1, 1, 2, 4, 8 and 16 times 
230 /^Jy/beam. The restoring beam is 4.4x1.7 mas, with major 
axis PA= -2.6° E of N. 



Analysis Guid^f] and the AIPS CookboolQ, respectively. 
After the initial amplitude calibration, parallactic angle 
correction (EVN) and Earth Orientation Parameter cor- 
rection (VLBA), the data were fringe-fitted in AIPS. At 
the first three epochs, when MAXI J1659— 152 was still 
brighter than 1 mjy, we also looked for polarized emis- 
sion in the VLBI data. We performed cross-hand fringe- 
fitting on a bright calibrator to remove the inter-channel 
R-L delay and phase offsets at the reference antenna. 
For the EVN data we used OQ208, whereas for the first 
epoch VLBA dataset we used NRAO530 . We then im- 
aged and self-calibrated OQ208 in Difmap (|Shepherd et al.l 
199|), and the self-calibrated OQ208 data were used to de- 
termine the antenna polarization leakage terms with LP- 
CAL in AIPS (for both the EVN and VLBA data). The 
phases of both the science target and the secondary calibra- 
tor were phase-referenced to J1707— 1415. The secondary 
calibrator, J1700— 1456, was imaged and self-calibrated in 
Difmap. This source has resolved structure on a scale 
of a few mas (see Fig. [TJ. In the higher resolution 
VLBA data, the average position of the maximum bright- 
ness peak in the brighter, South- Western component is 
a(J2000) = 17:00:03. 333340 and 8 (J2000)=-14:56:21. 97584. 
The statistical error from the four measurements is 170 /ias 



3 DATA REDUCTION AND RESULTS 



The EVN and the VLBA data were reduced according 
to the standard procedures outlined in the EVN Data 



3 http:/ /www. evlbi.org/user_guide/guide/userguide. html 

4 http:/ /www. aips.nrao.edu/CookHTML/CookBook. html 
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Figure 2. Naturally weighted 5 GHz EVN and VLBA images of MAXI J1659-152. The coordinates of the map centres are 
a(J2000)=16:59:01. 676959, 5 (J2000)=-15:15:28. 73200. The peak brightnesses are (in time order) 6.2, 3.9, 2.5, 0.9, 0.4 and 0.2 (non- 
detection) mjy/beam. The contour levels are powers of two times the 3a noise level, typically 170 /ijy/beam. The cross indicates the 
average VLBA position. Map parameters are listed in Table [l] 
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Table 1. Map parameters for Fig. [2] as well as modclfit results and brightness temperatures for all observations at 5 GHz. The columns 
are: observing array, observation start date [dd.mm.yyyy], Modified Julian Date of the midpoint of the observations, peak brightness 
[mjy/beam], noise [mjy/beam], restoring beam major and minor axes [mas], position angle [degrees] (note the last epoch map is a 
dirty map, therefore a restoring beam is not given), model-fitted circular Gaussian component flux density [mjy], radius [mas], and the 
corresponding (lower limit to the) brightness temperature. 
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1.1 


12.5 


VLBA 


06.10.2010 


55474.0 


0.87 
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-1.9 


1.07 
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VLBA 


14.10.2010 
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0.45 


< 0.7 


> 4.7 


VLBA 


19.10.2010 


55489.0 


<0.30 


0.057 


N/A 


N/A 


N/A 


<0.30 


N/A 


N/A 



in right ascension and 250 pas in declination. Accounting for 
the positional error of the primary phase-reference source, 
this increases to 320 /ias and 560 /ias, respectively. 

The image was imported back into AIPS where phase 
solutions were derived with CALIB using a solution interval 
of 30 minutes, and applied to the target. This step helped to 
improve the image fidelity of the target, as well as the astro- 
metric precision of the results. Note that we used the highest 
quality VLBA map to self-calibrate all epochs, except for 
the first epoch EVN data because the accurate coordinates 
of J1700— 1456 were not known during the correlation of the 
first experiment, therefore exact alignment of the VLBA and 
the first epoch EVN phase-centers is not possible below the 
~mas level. Choosing a 30 minute solution interval ensured 
that we corrected for a smoothly changing residual phase re- 
sulting from the imperfect troposphere model, rather than 
short timescale phase fluctuations. In the case of the EVN, 
we did not solve for Cambridge and Hartebeesthoek because 
they had low signal-to-noise ratios on J1700— 1456. In addi- 
tion, at the first epoch the secondary calibrator phase errors 
were higher than expected for the WSRT, likely because the 
a priori position had an error of about an arcsecond, which 
is a significant fraction of the phased-array WSRT beam. In 
that case the WSRT phases (only) were self-calibrated on 
the target, which had a comparable total flux density to the 
secondary calibrator at the first epoch. This additional phase 
self-calibration, using a 30 minute solution interval, did not 
significantly affect the final result. The final images of the 
target were made in Difmap, with no further self-calibration 
at the other epochs. 

In the following discussion we utilize the fact that MAXI 
J1659-152 had a flat spectrum during the first two EVN 
epochs. As part of a broadband monitoring campaign, the 
source was observed with the Karl G. Jansky Very Large 
Array (VLA) at multiple frequencies on 2010 September 29 
and October 1. On the former date the flux density was 
9.88 ± 0.30 mjy at 4.9 GHz and 10.03 ± 0.31 mjy at 8.5 
GHz, while on the latter date those were 10.29 ± 0.32 mjy 
and 9.74 ±0.30 mjy, respectively. The resulting spectral in- 
dices are 0.03 ± 0.11 and —0.10 ± 0.11, consistent with an 
index of 0. We note that there are no significant detections 
of polarization at those epochs, with upper limits of a few 



percent on the degree of polarization. Further details on the 
VLA data analysis can be found in van der Horst et al. (in 
preparation). 

3.1 Results 

The resulting VLBI maps are shown in Fig. O and the map 
parameters and modelfit results are listed in Table [T] At 
the first epoch on 2010 September 30 the source is well re- 
solved, as also clearly evidenced by the decreasing visibility 
amplitude with increasing baseline length seen in the uv- 
data, especially on the most sensitive baselines that include 
either Effelsberg or Westerbork. Fitting a single elliptical 
Gaussian model component to the uv-AaXa, gives a total 
flux density of 8.7 mjy and a characteristic size of 8 mas 
(well exceeding the uniformly weighted beamsize including 
all baselines), elongated at PA~132 degrees (measured from 
North through East). Analyzing the Westerbork synthesis 
array data taken during the VLBI observations, the total 
flux density of MAXI J1659-152 was 9.8 mjy, i.e. the single 
component contains most of the flux density from the target. 
To recover all of the WSRT total flux density, three model 
components were fitted to the data (see Table [T]). Note that 
the naturally weighted EVN images shown in Fig. [2] do not 
include the baselines to Hartebeesthoek, for better recon- 
struction of the extended emission. The uniformly weighted 
image including Hartebeesthoek data suggests there is ex- 
tension to the South-East on mas scales, in agreement with 
the modelfit results. This image contains a total cleaned flux 
density of 6.3 mjy. The naturally weighted image on the 
other hand has a total cleaned flux density of 9.6±0.1 mjy, 
recovering practically all of the integrated flux density. It 
shows an extension roughly to the East up to 20 mas. It is 
hard to judge from these data alone whether the apparent 
extension to the West-North- West is real or not, but this 
feature agrees in orientation with the major axis PA of the 
fitted elliptical Gaussian component. 

The higher resolution VLBA image on October 2 shows 
a partially resolved central component containing 6.3 mjy. 
Fitting a single elliptical Gaussian component resulted in a 
major axis of PA~115 degrees, close to the single-component 
fit EVN value. There is extended emission on both sides of 
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the source roughly in this direction. This extended emis- 
sion is only detected on the shortest spacings. But since it 
does not disappear even with further point source model 
phase self-calibration, and the overall structure is similar 
to the EVN result, we conclude that the extended emis- 
sion is likely real. The source was therefore fitted with three 
circular Gaussian components, with a total flux density of 
8.8 mjy. By the second EVN observation, on October 4, 
the total flux density of MAXI J1659— 152 had decreased to 
2.95 mjy according to our WSRT measurement. We fitted 
a circular Gaussian brightness distribution model directly 
to the uu-data. The entire integrated flux density was re- 
covered this way in a single component (consistent with the 
WSRT value within the errors), with a radius of 1.1 mas. 
MAXI J1659-152 was detected in two more VLBA obser- 
vations. Although the source looks compact in the images 
from epochs 3 and 4, the source was still resolved according 
to our modelfit results (see Table [l}. At the fifth epoch the 
fitted Gaussian component size was smaller than the the- 
oretical resolution limit of our array, therefore we consider 
the minimum resolvable angular size as an upper limit (see 
lLobanovll2005l ). MAXI J1659-152 faded below the detection 
level of 300 /ijy/beam (5a) on 2010 October 19. 

No polarized emission was detected at any of these 
epochs. The first EVN epoch, with the highest peak in total 
intensity and lowest noise in polarized intensity, provides the 
strongest constraint. The 150 /xjy/beam upper limit (5a) 
in polarized intensity corresponds to a 2.5% upper limit in 
fractional polarization. 

Using the higher resolution VLBA 
data, the derived average position for 
MAXI J1659-152 is a (J2000)=16:59:01.676891 and 
S (J2000)=-15:15:28. 73237. The error in declination 
is 220 pas, consistent with the measured errors in the 
secondary reference source position. In right ascension 
the error is much larger, 690 /xas. The peak in the (lower 
resolution) first epoch EVN map is about 2 mas East of the 
average position. This is unusual since due to calibration 
errors at low elevations, one expects to see larger errors i n 
declination than in right ascension (cf. lYang et al.l 12012). 
Therefore, this is likely the sign of a source structural effect. 
The October 4 EVN position is within the 2a error. The 
total errors in a and S, including the uncertainty in the 
primary phase-reference source position, are 740 pas and 
550 //as, respectively. 



4 DISCUSSION 

4.1 A compact jet in MAXI J1659-152? 

The images in Fig. [2] clearly show resolved or partially re- 
solved structures at the first three epochs, indicative of an 
outflow. The maximum brightness temperatures resulting 
from fitting circular Gaussian components to the ww-data at 
each epoch range from ~ 10 7 — 10 8 K (see Table [T]). These 
high brightness temperatures, together with the observed 
flat radio spectrum (see Sect. 3), are indicative of optically 
thick synchrotron radiation, which is also consistent with the 
observed low fractional polarization (<2.5%). An optically 
thick synchrotron spectrum may originate in partially self- 
absorbed compact jets (for details see Appendix). In fact, 
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Figure 3. The jet intensity profile with distance from the cen- 
tral engine. Both dS/dr and r are in arbitrary units. r core is the 
distance between the jet origin (~the compact object) and the 
brightest peak of the compact jet. It demonstrates that the ob- 
served 'radio cores' of microquasars and AGN do not exactly co- 
incide with the location of the black hole. d COTe is the FWHM size 
of the compact jet profile (approximately the size of the source 
obtained by fitting a Gaussian profile to the uv-data). The model 
is scaled so that the optical depth is unity at distance r=20, at 
reference frequency vq. Jet parameters are m=l, n=2, s=2 (k r =l, 
fc T =3). For an explanation of the different jet parameters, see the 
Appendix. 



such partially self-absorbed compact jets have been resolved 
with VLBI in the BHXRB candidate SS433 (|Paragi et all 
Il999h. and the confir med BHXRB sy stems GRS 1915+ 10 5 
ijDhawan et all l2000h and Cyg X-l l|Stirling et al l l200ll ). 
Various other observational evidence points to the fact that 
in general, BHX RBs produce compact jet s in the canonical 
hard X-ray state ^Fender et al.ll2004l2009l ). Our highest res- 
olution data do not show a resolved core-jet structure, but 
the images from the first two epochs clearly show extensions 
roughly to the East and to the North- West, that may be re- 
lated to an outflow. The model-fitting results are also consis- 
tent with an elongated structure in a similar direction. One 
explanation for the overall fuzzy appearance of the 'core' 
could be angular broadening of the images due to inter- 
stellar scattering. Definitive evidence for scatter-broadening 
would be multi-frequency measurements that showed & v~ 2 
dependence of the observed source size. However, the nearby 
secondary calibrator source J1700— 1456 showed no evidence 
for scatter-broadening. 

We have two other pieces of circumstantial evidence for 
supporting the self-absorbed jet nature for the bulk of the 
radio emission. Curiously, MAXI J 1659— 152 appeared to be 
the most compact at the last detection on 2010 October 14, 
when the total flux density was less than 1 mjy. At ear- 
lier epochs (and especially with the larger EVN beam) a 
larger overall source size was measured because of the ex- 
tended emission (outflow). The radius of maximum bright- 
ness temperature, and consequently the measured size of 
the compact jet region, d epends on the radio luminosity 
|Blandford fc KonigifnmSI) . The measured size values from 
Table [TJ are plotted in the left panel of Fig. [5] as a function 
of flux density. While the scatter is significant, there seems 
to be a trend of increasing angular size with increasing flux 
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density, as expected for compact jets. The different array 
configurations of the EVN and VLBA observations proba- 
bly affect the size measurements, so we cannot compare the 
size-flux dependence with model predictions. The important 
point here is that changes in total jet power or bulk Lorentz 
factor will change the source flux as well as apparent size: 
initially the source brightens and grows, but as the flare 
fades it shrinks - this is not expected in cases where the ra- 
dio emission is related to discrete ejecta or a shell of matter 
that was ejected from the system. 

A strong dependence of size with flux dens ity was previ- 
ously found in another microquasar, Cyg X-3 (| Newell et al.l 
If 998h . In that case the authors claimed to have observed 
superluminal expansion and contraction on timescales of 
hours. Another similarity is that while a collimated ra- 
dio jet is clearly present during outbursts of Cyg X-3 (e.g. 
iMioduszewski et alj|2001f ). the optically thick base morpho- 
logically does not resemble a core-jet structure. Instead, 
it looks like the 'fuzzy' core of MAXI Jf 659- 152. In the 
case of Cyg X-3, the scatter-broadening effect is well estab- 
lished, with a scattering d isk size of ~18 mas at 5 GHz 
l|Mioduszewski et al.l l200lh - However, the explanation of 
scattering is inconsistent with the small size we measured 
in the last epoch. This same argument (i.e. the measured 
source size being smaller than t he scattering size) could also 
imply that lNewell et ail (119981 ^ may not have measured ac- 
tual physical sizes, thereby weakening the claim for appar- 
ent su perluminal expansion in Cyg X-3 (|Mioduszewski et al.l 
l200lf ). In the case of MAXI JI659-I52, a scattering screen, 
if present, must be well localized at the target position be- 
cause scatter-broadening is not observed in the nearby cali- 
brators. 

Besides the apparent correlation between the source 
core flux density and its size, we compare the dependence 
of the apparent core position with the variable core size, be- 
cause in compact jets these are intimately related. In Fig. [3] 
we show the intensity profile for a simple compact jet model 
(see Appendix). The observed radio core position as well as 
its apparent size change with observing frequency. This is 
because at different frequencies the optical depth is differ- 
ent; a n effect known a s the frequency dependent core-shift in 
AGN <|Lobanovlll998h . The true position of the black hole i s 
not easy to determine with high accuracy. lHada et al.l (1201 lh 
used mult i- frequency observations to pinpoint the location 
of the jet base in M87, and they concluded that the su- 
permassive black hole is within 14-23 Schwarzschild radii 
of the 43 GHz radio core (a projected angular separation 
of 41 fias). The only similar measurement for Ga lactic mi- 
croquasars was made in SS433. |Paragi et al.l (|l999h used the 
well-established kinematic model of the ballistically moving, 
optically thin jet components to locate the binary system in 
between the approaching and receding sides of the optically 
thick radio core components. The separation between the bi- 
nary system and the observed radio core position decreased 
with increasing frequency, i.e. they observed the frequency 
dependent core-shift in a Galactic microquasar. 

Because we observed MAXI J1659— 152 at a single fre- 
quency, we cannot probe the frequency dependent core shift. 
However, any change of optical depth during an outburst 
(due to e.g. an increase in relativistic particle number den- 
sity and/or an increase in the magnetic field strength) will 
result in a similar effect when we compare core positions and 



sizes at various epochs at a single frequency. This will result 
in a larger separation of the observed peak radio emission 
and the compact object, and at the same time a larger ob- 
served size for the compact jet. This effect can be probed 
with VLBf at a single frequency, even if, as mentioned above, 
the core shift, r core , cannot be measured in most cases. By 
simply scaling our model with the initial optical depth, we 
find that r coro is roughly proportional to the FWHM of the 
jet profile, d C orc, in which case the difference of peak posi- 
tions Ar will be proportional to the difference of the mea- 
sured core sizes Ad between the various epochs. In Fig. [4] we 
plot the relative peak positions of MAXI J1659— 152. Note 
the apparent positions in right ascension continuously drift 
from roughly East towards the West at the various epochs, 
indicating either proper motion or varying source structure. 
Because the first epoch map clearly shows an extension to 
the East, it is plausible that we see a shrinking compact 
jet as the source flux density decreases. In the right panel of 
Fig.[5]we plot the difference in core position versus difference 
in core size, with respect to the first epoch. The data show 
a smooth change in source size along with a smooth change 
in core (more precisely, compact jet peak) position in a way 
that is compatible with the expected relation between r COIC 
and dcorc for compact jets. Monitoring VLBI observations 
of similar transients at multiple frequencies can potentially 
confirm these effects: the core size and position change with 
luminosity ( luminosity dependent core-shift), as well as with 
observing frequency (frequency dependent core shift, shown 
in Fig. O . 

On the other hand, one may interpret the positional 
change as source proper motion. Here we take a conser- 
vative approach. Because the first epoch was clearly af- 
fected by source structure and the time baseline is very 
short, we give a 2a upper limit of 1380 fias in right as- 
cension and 440 fias in declination; in twelve days (epochs 
2-5), these correspond to proper motion upper limits of 
fi a < 115 ^tas/day and (is < 37 /ias/daj|j. To convert 
these to velocities, we have to know the source distance, 
which is not strongly constrained, as is typically the case 
for X-ray transients. Based on the relation between opti- 
cal outb urst amplitude and the o r bital period in s imilar 
systems (|Shahbaz fc Kuulkers|[l998r ). lKuulkers et all (|2012h 
estimated a di stance of 7±3 kpc for M AXI J1659-152. As 
pointed out by iMiller- Jones et alj (|201ll ). this does not take 
into account the effect of orbital inclination to the line of 
sight; at high inclinations the observed brightness of the 
disc would be underestimated, and consequently, the de- 
rived distance would b e an overestimate. In a subsequent 
paper, iKuulkers etafl (|2013T ) derived an updated distance 
of 8.6±3.7 kpc, using the original amplitude-period rela- 
tion. Correcting for the inclination angle of the orbit, the 
distance would be lowered to 7.1±3.0 kpc at i — 65°, and 
4. Oil. 7 kpc for i = 80°. The corresponding heights above 

5 We note that in order to minimize the structural effects in right 
ascension, one could use the last two VLBA detections to con- 
strain the proper motion in right ascension, when the source was 
most compact. Although the time baseline is somewhat shorter 
(9 days), the measured statistical errors in right ascension for the 
secondary calibrator are also smaller. Using these we would arrive 
at a very similar proper motion constraint in right ascension as 
in declination. 



8 Z. 



Paragi et al. 



5470 



5475 



5480 



5485 



MJD-50000.0 



E 



5470 



5475 



5480 



5485 



MJD-50000.0 



Figure 4. Map peak position relative to the map centre in right ascension (a cos <5; left panel) and declination (<5; right panel), as 
measured at 5 GHz at the various epochs. The EVN measurements are indicated by triangles, the VLBA ones with squares. There 
is a clear, smooth change in right ascension during the VLBI monitoring programme, most likely related to structural changes in the 
emission. The increasing negative offset in right ascension shows a gradual shift of the source peak emission to the West. There is no 
clear variation in declination, and the measurements agree within the errors. 
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Figure 5. The left panel shows the source size of MAXI J1659-152 (circular Gaussian component fit; we used the dominant core 
components at the first two epochs) as a function of its flux density at 5 GHz (triangles: EVN; squares: VLBA). In the right panel 
we show the observed 'core' position change (A r) versus the measured source size change (A d) relative to the first epoch. From a 
simple scaling of the opacity in our compact jet model, it is expected that the position of the peak emission and the FWHM width of the 
compact jet scale linearly with one another. Therefore, by analogy with the frequency-dependent core shift (see Fig. [3]), one may measure 
a time-dependent core shift at a fixed frequency during an outburst, and compare it to the observed change in core size. The data seem 
to show the expected trend. Admittedly, the number of measurements is small, and the errors are large (triangle: EVN; squares: VLBA). 



the Galactic plane are z = 2.4 ± 1.1 kpc, z = 2.0 ± 0.9 kp c 
and z = 1.1 ± 0.5 kpc, respectively jKuulkers et al.|[2013h . 
Our VLBI images do not show a highly beamed source struc- 
ture, with a very compact core and one-sided ejection, and 
they do not show symmetric double-sided ejecta in the plane 
of the sky either. Assuming that the jets are perpendicular 
to the accretion disk, this would mean that the orbit of the 
binary is likely moderately inclined to the line of sight, pre- 
ferring the ~7 kpc value. Finally, by carefully taking into 
account various arguments, Ijonker et al.l |2012i 'l estimated a 
source distance of 6 ± 2 kpc, although this error bar may be 
underestimated given the number of assumptions involved. 

Assuming a distance of 7 kpc, our proper motion con- 
straints correspond to upper limits on the space veloc- 
ity of <1400 kms -1 in right ascension and <440 kms -1 



in declination. In the high inclination angle case (corre- 
sponding to a distance of 4 kpc), these numbers would 
change to ~800 kms -1 and -~250 kms -1 , respectively. 
These values are not particularly constraining, since there 
are no measured blac k hole kicks as large as 1400 kms -1 . 
lYamaoka ef~aH (|2012T ) argued that MAXI J1659— 152 is a 
runaway microquasar, just like XTE J1118+480, which has 
a velocit y of 145 kms -1 with r espect to the Lo c al Sta ndard 
of Rest |Mirabel et al.l 120011 ). iKuulkers et~al] (|2013l ) con- 
cluded that this was a likely scenario for MAXI J1659— 152 
because shorter period (less massive) systems should have 
the highest runaway velocities; this scenario would also ex- 
plain the large distance from the Galactic plane. Our twelve 
day time baseline was too short to significantly constrain 
the velocity. With the accurate coordinates measured during 
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Figure 6. RXTE PCU2 count rate (2-20 keV, left) and hardness-intensity diagram of MAXI J1659-152 during its outburst in 2010. 
The hardness is defined as the ratio of count rates measured in the 6-15 keV and the 2—6 keV energy ranges. The black bars indicate 
the times of the VLBI observations while the black circles show the RXTE measurements closest in time to the VLBI observations. 



this outburst, future VLBI observations in the next activ- 
ity cycle will give an accurate space velocity measurement 
for MAXI J1659— 152. For future X-ray transients with a 
more prolonged flaring activity, it will be possible to use 
this method to probe the runaway microquasar scenario if 
the transient is within a few kpc. Further improvement on 
the positional accuracy can be obtained if one can find a sec- 
ondary calibrator closer to the target, especially one within 
the primary beam of the telescopes. But as we have seen in 
here, varying source structure will likely be a limiting factor 
in most cases. 



4.2 X-ray properties and radio structure 

The observed radio structure of MAXI J1659— 152 is consis- 
tent with a compact jet, which is typical for the canonical 
hard state of BHXRBs. We will now consider the X-ray prop- 
erties of the transient. The RXTE lightcurve and HID (6—15 
keV/2— 6 keV) are shown in Fig. [U black points show mea- 
surements closest to the VLBI epochs. The source already 
left the canonical hard state and was in the hard intermedi- 
ate state (HIMS) by the time of the first EVN epoch. The 
presence of a compact jet in the HIMS state is not unprece- 
dented; Cyg X-l is known to show a compact jet in this state. 
During the evolution to the soft intermediate state (SIMS), 
radio qu enching was reported from M AXI J 1659— 152 on Oc- 
tober 8 (|van der Horst et al.ll2010b|) and tvpe-B QPOs ap- 
peared on October 12 (|Munoz-Darias et al.ll20Tll ), which are 
usually followed by relativistic jet ejections. Note that radio 
quenching and subsequent bright, discrete ejecta are typical 
of BHXRBs during their transition to the soft state. How- 
ever, while our measurements confirm the radio quenching, 
we do not see discrete ejecta in MAXI J1659— 152. Instead, 
the source simply faded below our detection limit. Because 
of the relatively rapid cadence (~2 days) of our observations 
during the radio bright state of the source, and because no 



other flaring events were reported from total flux density 
monitoring campaigns, it is quite unlikely that we missed an 
ejection event. We note also that our final VLBA observa- 
tion was taken when the source was in its softest X-ray state 
(see Fig. HU, and by that time the RXTE count rate was well 
into its decay phase, therefore we do not expect significant 
radio flaring after this time. The lack of ejecta may indi- 
cate that MAXI J1659— 15 2 underwent a failed trans ition. 
However, as was shown by iMunoz-Darias et al.l l|20Tll ). the 
X-ray state transition was fully completed, according to the 
well established criteria: the total rms variability went down 
to almost 1%, unlike H1743— 322. It is true however that 
the soft state in MAXI J1659— 152 was very short-lived and 
appeared significant l y hard er than that in most BHXRBs. 
IMunoz-Darias et al.l l|201ll ) proposed that this may be due 
to the high inclination of the system, in which case we see 
less of the soft disc emission while the emission from the 
spherical, non-thermal corona should be largely inclination 
independent. A similar lack of bright, transient ejecta was 
reported fo r the transition of C yg X-l to the soft state in 
2010 June (jRushton et al.ll2012l ). However, Cyg X-l is also 
peculiar in that it has a non-canonical soft state (signifi- 
cantly higher fractional rms variability than typical tran- 
sient BHXRB outbursts). VLBI observations of many more 
transient events are necessary to draw general conclusions 
about the formation of transient ejecta, whether they appear 
in all soft-state transitions, and where exactly the 'jet-line' is 
located in the hardness-intensity diagram (which may vary 
from system to system , and even within a given system; 
iMiller- Jones et alllioT^ V 

Ijonker et alj (|2012h disc uss the well-known radio-X-ray 
corre l ation in BHXRBs (e.g. ICorbel et al"1l2003l : iGallo et all 
120031; IJonker et al.ll2004l ). During the canonical hard state, 
the observed correlation is Lr oc L^ 6 ; more recently it has 
been shown that a number o f systems follow a much steeper 
correlation with Lr oc L^ 4 l|Coriat et al.l [20T1I ; IGallo et al.1 
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12012ft . At the beginning of its o utburst, while the lu minos- 
ity was in excess of 10 erg s" ' Mjonker et~ai] (|2012h found 
MAXI J 1659- 152 to follow the steeper correlation, but to 
be significantly more radio luminous at a given X-ray lumi- 
nosity than the rest of the sources on the steeper correla- 
tion track. To explain this discrepancy, they proposed that 
the source either falls on the outlier branch, but with sup- 
pressed X-ray luminosity owing to the high inclination, or 
that the radio emission was possibly dominated by optically 
thin ejecta at early times, incre asing the average rad io lumi- 
nosity. In contrast, in figure 9 of lCorbel et all (|2013h . MAXI 
J1659-152 does not appear to be significantly radio overlu- 
minous for the steep er correlation track. A re-examination of 
the original data of IJonker et al.l ((2012) revealed a discrep- 
ancy in the assumed distance of H1743-322, the single source 
that best defines the steeper correlation track (|Coriat et al.l 
l201ll ). Assuming instead the accepted distance of 8 kpc for 
H1743-322, MAXI J1659-152 would no longer appear signif- 
icantly radio overluminous (P. Jonker, priv. comm.). How- 
ever, this does not address an additional motivation for the 
possibility of optically-thin radio flares namely the reported 
detection of a high level of line ar polarisation from a prelim - 
inary analysis of WSRT data (|van der Horst etafll2010al ). 
Reanalysis of these WSRT data however showed no signif- 
icant linear polarization (see Sect. 3), in agreement with 
our VLBI results, which show no evidence for optically thin 
ejecta. The small observed physical size, the flat spectrum, 
the lack of significant proper motion and the low linear po- 
larization suggest that the radio emission was dominated by 
a compact jet. 



5 CONCLUSIONS 

We presented a series of VLBI (realtime e-VLBI at the EVN 
plus VLBA) monitoring results for the BHXRB candidate 
MAXI J1659-152. The VLBI maps show evidence for an 
outflow on a scale of 10-mas. The compact jet scenario fits 
the data well: the low polarization, flat spectrum as well as 
the observed source size and luminosity dependent core shift 
all support this interpretation. During the full phase transi- 
tion from the HIMS through the SIMS to the soft state, we 
do not see evidence for the luminous, discrete ejecta that are 
typical for BHXRBs in this state, possibly indicating that 
the outflow rate was not sufficie nt to form powerful sh ocks, 
as seen in the case of Cyg X-l. As lYamaoka et al l \201$ ) and 
iKuulkers et~all (|2013l ) pointed out, MAXI J1659-152 could 
potentially be a runaway microquasar. Further VLBI obser- 
vations during another outburst will be necessary to measure 
the proper motion of the system with high accuracy. With 
the method presented here, it should be possible to pro- 
vide useful constraints on black hole transient proper mo- 
tions even during a single outburst if the outburst lasts long 
enough. The runaway microquasar hypothesis could then be 
tested. 
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APPENDIX A: COMPACT JET MODEL 

A conical jet model to describe the intensity profile of the 
radio jets in the ga l actic XRB SS433 was developed by 
iHiellming fc Jonstonl i|l98Sl ). hereafter HJ88. Here we follow 
HJ88, but instead of the large-scale jet we will concentrate 
on the properties of the optically thick inner core, also taking 
into account more general considerations for the variation 
of magnetic field with distance, as done for case of AGN jet 
models (see below). Let us assume a conical and freely ex- 
panding jet with radius proportional to the distance r to the 
central engine. Synchrotron radiation is emitted by electrons 
moving in a tangled magnetic field with average strength B. 
We assume that the energy spectrum of the radiating elec- 
trons - filling the entire volume of the jet - has a power-law 
form, iV(7 e ) ~ 7,T S , where s is the energy spectral index, 
which is generally thought to be ~2 (1.0 < s < 3.0 will be 
considered in the present paper). The optically thin spec- 
trum of the so urce will have a spect r al ind ex a — (1 — s)/2 
(S et e.g. iRvbicki fc Lightman! <|l979h >. The magnetic 
field and the particle density are assumed to depend on r 
as B oc r~ m and N oc r~ n ; for a conical geom etry, m can 
vary between 1 and 2 l|Blandford fc Ree3ll974h . If m = 1 
and n = 2, the energy density of the magnetic field and the 
kinetic energy density of the electrons are in equi partition at 
any r along the flow (|Blandford fc Konigllll979h . All quan- 
tities determined below are expressed in the emitting frame 
of the fluid. 

The optical depth to syn chrotron radiati on through the 
jet axis can be written as (cf. iLobanovll 19981 ) 

T(r) = T () (r/ro)- M£+1) Wo)- (£+1) , (Al) 

where e = (s + 2)/2 and k r = (em + n - l)/(e + 1). For 
the equipartition values of TO = 1 and n — 2, k r = 1 
independent of the energy spectral index. In this case the 
optical depth scales as r - ' 6 "^ 1 ' and the location of the ra- 
dio core is inversely propor tional to the observing frequency 
l|Blandford fc Ko nigl 1979). The actual value of the constant 
ro depends on the selection of ro and 6, the jet inclination 
to the line of sight (assumed to be constant). 

Following HJ88, the contribution to the flux density of 
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a portion of the jet of length dr can be written 

^ = Csintf {y/vof 2 (r/r ) (m+2)/2 (1 - 0€«*(t), 
ar 

(A2) 

where £ con = a + br — ct 2 ... is the geometrical correction 
function for conical jets and a,b, c . . . constants are given by 
iHiellming fc Hani (| 1995T ) . The actual value of the constant C 
is determined by the opening angle of the flow, the magnetic 
field strength and number density of electrons, etc., at a 
reference distance ro- Note that Eq. (IA2|) and Eq. (|A1|I are 
more general than that given originally by HJ88, since m 
and n are not fixed. These values, together with the energy 
spectral index s - through k r - determine the shape of the 
jet intensity profile (see Fig. [3]). 

The jet intensity first rises sharply with r and reaches 
its maximum at r = r corc (at a given frequency). Within 
a few r corc its flux density will decrease significantly. One 
can also determine the brightness variation in the core with 
frequency. The optically thi ck limit of th e integrated flux 
density from Eq. |(A2j gives l |Konidll98lh 

Qcorc = 5/2 - (m + 4)/(2fe r ) (A3) 

While the peak of the jet intensity profile is roughly pro- 
portional to the frequency, the integrated spectrum of the 
core remains flat or slightly inverted, as typically observed 
in radio- loud AGN. The spectral index can vary between 
—0.5 < a < 1.5, with the equipartition value of a — 0.0 . 

In our model k r is assumed to be constant throughout 
the region of interest (in the core at a given frequency). If 
k r changes smoothly with decreasing r (consequently, at in- 
creasing v), then Eq. (|A2|) is approximately valid only in a 
frequency range Au — V2 — vi, for which k r is determined. 
In the inner jets (at v 3> 1 GHz) however, the geometry of 
the flow may be dif ferent from conical due to external pres- 
sure gradients (e.g. iGeorganopoulos &i Marscherill996l ) and 
k r may change very fast. Clearly, our simple model does 
not describe this scenario well. External pressure gradients 
may be caused for example by a dense, optically thick stel- 
lar wind. As most of the jet environment is expected to be 
ionized, this results in additional free-free absorption, and 
the jet intensity profile will depend also on the distribution 
of the absorbing medium. 

This paper has been typeset from a TeX/ IMgX file prepared 
by the author. 



